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An energy-dependent and a set of single-energy partial-wave analyses of 
Qv . Trd elastic scattering data have been completed. Amplitudes are presented 

for pion laboratory kinetic energies up to 500 MeV. These results are com- 
'^ . pared with those found in other recent analyses. We comment on the present 
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^^ . database and make suggestions for future work. 
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I. INTRODUCTION 

The vrd system is an important component of the more general ttNN problem. At 
intermediate energies, the ird and NA channels provide much of the inelasticity associated 
with pp scattering |jl]]. We have previously analyzed the pp — > pp and vrrf — ;> pp 
reactions. Here we give the results of our analyses of ird elastic scattering data. 

A number of calculations have resulted in a reasonable qualitative (and over some regions 
quantitative) description of red scattering data |JI|. The more theoretical approaches have 
had difficulty in describing all observables |^, and have generally concentrated on limited 
kinematic regions. In addition, several partial-wave @]- P] analyses have been performed. 
These analyses have found some motivation from claims [§ of dibaryon resonances in pp 
and nd —>■ pp scattering reactions. Comparisons with data have been complicated by the 
occasional appearance and subsequent disappearance of sharp structures in some polariza- 
tion observables. We plan to examine these questions through a coupled- channel analysis 
of pp and ird elastic scattering combined with ird -^ pp. However, as a first step, we have 
analyzed nd elastic scattering separately. 

The present analyses include measurements to 500 MeV in the laboratory kinetic en- 
ergy of the pion. (This corresponds to laboratory kinetic energies between 287.5 MeV and 
1287 MeV in the pp system.) The value of -^i varies from 2.015 GeV to 2.437 GeV, and 
spans the range of energies typically associated with dibaryon candidates. This reaction 
displays the same "resonancelike" behavior found in partial-wave analyses of pp — > pp and 
ird — > pp scattering data. The interpretation of this behavior generally arouses strong reac- 
tions, either for or against the existence of intermediate dibaryon resonances. In the present 
work, our goal is simply a refined understanding of the nd elastic scattering amplitudes. 
Interpretations will depend upon the results of our larger joint analysis of the pp and nd 
channels. 

In the next section, we will make some comments on the database used in this analysis. 
In Section HI, we will outline the general formalism for vrd elastic scattering. Methods 



used in the partial-wave analysis will be discussed in Section IV. Our main results will be 
presented in Section V. Here we will also compare with the available data and other recent 
analyses. Finally, in Section VI, we will summarize our findings and make some suggestions 
for future investigations. 

II. THE DATABASE 

Experimental studies of vrd elastic scattering began to produce results in the 1950s. At 
this time, the first measurements of total and differential cross sections became available. 
The trend of vrd elastic scattering data accumulation since 1952 is displayed in Fig. 1. The 
rapid increase in the number and type of measurements in the early 1980s was motivated 
by a growing interest in the problem of exotics. This reaction was expected to give further 
information on the existence (or non-existence) of dibaryon states suggested in analyses of 
A^A^ elastic scattering data. Numerous high-quality deuteron polarization measurements 
were made. The total database more than doubled, and several partial-wave analyses ^- 
P] were carried out at this time. The present study has utilized a set of data which is 
significantly larger, and covers a broader energy interval. 



Our total set of experimental data |I0|^ [131 (1362 points) includes measurements of 
the TT+d (57) and vr^d (67) total cross sections (cr^), Tr+d (516) and 7r~d (236) differential 
cross sections (da/dQ) with unpolarized targets, the deuteron vector analyzing-power (iTn) 
for TT+d (280) and Tr^d (5), the deuteron tensor analyzing-power (T20) for 7r"^d (42), the 
combined deuteron tensor analyzing-powers (r2i and T22) for Tr^d (123), and the laboratory 
deuteron tensor analyzing-power (tgo'') for 7r+d (30). We have also included 6 unpolarized 
total elastic cross sections {<Jt)- Energy-angle distributions for da/dVt, iTu and T20 are 
given in Fig. 2. Most of the vr+d data are concentrated at low energies; the 7r~d data tend 
to span a wider energy range. 

As shown in Table I, we have removed ar (2), da/dft (88), a^ (3), iTu (6) and t2o'' 
(35) data corresponding to 9% of the total. These measurements were the source of serious 



conflicts within the database and were not included in the analysis. For instance, many 
of the tgo^ measurements (33 points) were produced by the Ziirich group [0. These were 



found to be in conflict with a number of independent measurements at LAMPF p7|,p8|,|56 



TRIUMF [|g,g and PSI |3 



III. FORMALISM 

The relations between partial-wave amplitudes and observables have been given in a 
number of previous theoretical and phenomenological studies. Many of the these relations are 
given in the work of Grein and Locher ||58|. They have been included here for completeness 
and to deflne our notation. Due to parity conservation, there are 4 independent helicity 
amplitudes for this reaction. Thus, for a reconstruction of the scattering amplitude at flxed 
values of the energy and scattering angle, one requires 7 independent measurements. The 
amplitude, HaplO), is labeled by the deuteron helicities {a and (3) in initial and flnal states. 
Here the angle 6 is the center-of-mass scattering angle of the outgoing pion. Our notation 
for the helicity amplitudes |Q is given below: 
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H+^^H, = -Y: [{J+l) Tl, j_, + J Tj^, j^, + (2 J+ 1) ry , + 2 V^J (J+ 1) Tj_, j^,] d(„ 
^ j>i 



-f^-i-n = 



^2 = -^ E [\/2(J + l) (Tj^, j^, - Tj_, j_,) + V2 Tj_, 



J+lJ "10' 



:i) 



j>i 



H^_^H, = -Y: [(J+l) Tj_, j_, + J Tj^, j^, _ (2 J+ 1) T/, + 2 v^J (J+ 1) Tj_, ,^J d'U 
^ j>i 



Hoo ^H,= J2[J TU j_, + ( J + 1) T/+1 j^, - 2^. J (J + 1) TU j^,] rf, 

J>0 



J 

00; 



where the d^^ are reduced rotation matrices. The following symmetry relations 



Ha 13 
Ho 13 



-Hpo 



(2) 



are also obeyed by the above helicity amplitudes. The partial- wave amplitudes, T^^, ^^ , are 



labeled by the values of L^ corresponding to the vrrf flnal state, and L'^ for the -Kd initial 



state. In the next section, and in our figures, we use the notation ^L}, and ej for the 
amphtudes ^(J - l)j - 3(j + i)^ as in Ref. §. 

The various observables for vrd elastic scattering are given in terms of hehcity amphtudes 
5|,|59| below: 



Io=C = 2 \H^\' + 4 \H2r + 2 \H,\' + \H^\\ (3) 

da 

rSf = Anar, I In sin 9 dO 





aT = Anag [2ImHi{0) + ImH^iO)], 

zTn = ^t?? = -VQ Im[H;{Hi - H^ + H^)]/Io, 
T20 =C = V2 {\H,f - iH^f + m' - \H,f)/Io, 

T21 = -t"? = -VQ Re[H;{H, -H,- H,)]/Io, 
T22 = 4l = V3 [2 Re{HlH,) - {H^fj/Io, 
tl', = 3{\H,f-\H,f)/Io, 



t\l = 3Re[H*{H, + H,)]/I, 
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ztll = 3 V2 Im{HlH,)]/Io, 
tl\ = 3[Re{H*,H^) + \H2f]/Io, 
t\\, = 3 [-Re{H;H^) + \H2f]/Io, 
till = -V3 Im[H;{Hi -H,-2 H^)]/Io, 



it^i = -3 Im{H*Hi)/I^ 



0, 



itl\^ = -3 Im{H;H4)/Io, 

til = 3^2 Im{H;Hi)/Io, 
4^2 = -3 72 Im{H;H3)/Io, 



^20 = (l-f^ll 



4 \H2\^+\H^\^-2 |i/4|')//o, 



tf^ = VI Re[H;{Hr -H^ + 2 H,)]/Io, 
tl° = V6[Re{H*Hs) + \H,f]/Io, 
til = 3 [Re{HlH^) - lifslV^o, 

t^i = 3 72 Re{H;Hi)/Io, 



,21 
''2-2 



-3 V2 Re{H;H^)/k 



0, 



^22 — 3 l-f^ll //o, 
^2-2 = 3 \H-i\ /Iq. 



In these relations the factor ag is equal to | (^) , where k.,^ is pion momentum in 
the center-of-mass frame. Here the t's are spherical harmonics ||6^. The superscript * 



denotes complex conjugation. We have followed the Madison convention [^ in defining 
a coordinate system. The observables da/dfl, a^, and ar are the usual unpolarized cross 
sections, a'^^^"'^ gives the spin-aligned total cross section. The other observables involve 
either the polarization of the initial or final deuteron (iTn, T20, T21, and T22) or both. The 
following relationships are valid due to parity conservation and time reversal invar iance: 
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We also require some combinations of deuteron tensor analyzing powers measured by 



experimentalists. In Ref. |^5[, T21 and T22 are defined by 

[1 



T22 



T20 + ^22 , 



(5) 



1 _ 1 /I T. T. It. 

'?"2i — T21 + - T22 — - y - T20 + T21 + - T22. 



The laboratory deuteron tensor analyzing-powers expressed H3] in terms of CM variables 



are 



4o' = ^"'"''^^''"^ T20 + 2 ^^ sin ^« cos 9n T,, + ^ sin^ 9n T,,, (6) 

4f = J I sin 0R cos ^,j T20 + (1 - 2 cos^ ^^) T21 - sin Or cos ^,j T22, 

4f = ^/fsi^'^^ ^20 - sin^^cos^^ n, + i±^^ r22, 
where On is the deuteron recoil angle in the lab frame. 



IV. PARTIAL- WAVE ANALYSIS 

The energy-dependence of our global fit was obtained through a coupled-channel K- 
matrix form, in order to ensure that unitarity would not be violated. The "inelastic" channel 
was non-specific but included to account for the coupled pp and A^A reactions. For single- 
channel nd states (for example, ^02) this resulted in a 2x2 matrix, and for the spin-coupled 
states (for example, ^P2, €2, ^-^2) we had a 3x3 matrix. The matrix elements were expanded 
as polynomials in the pion laboratory energy, and an appropriate phase-space factor was 
included in the ird elastic elements to ensure proper threshold behavior. This analysis 
included 21 searched partial- waves and 66 varied parameters. Amplitudes with J < 5 were 
considered. The solution gave a x^ of 2743 for the 1362 data and 333 experiments below 
500 MeV. Coulomb modifications of the phase-space factors were attempted but discarded 
in the final fit; the sensitivity to such refinements appeared minimal. 

Single-energy solutions were produced up to 300 MeV, using a binning width of 
10±5 MeV. We used mainly the energy values chosen in previous single-energy PWAs [Q- [0. 
Starting values for the partial-wave amplitudes, as well as their (fixed) energy derivatives, 
were obtained from the energy-dependent fit. The scattering database was supplemented 
with a constraint on each varied amplitude. Constraint errors were taken to be 0.02. This 
was added, in quadrature, to 5% of the amplitude. Such constraints were essential to prevent 
the solutions from "running away" when the bin was relatively empty of scattering data. 
These errors were generous enough that they afforded little constraint for those solutions 
where sufficient data existed within the bin. 



Single-energy analyses are done in order to reveal "structure" which may be missing from 
the energy-dependent fit. Little compelling evidence was found for such structure. Results 
of the single-energy analyses are summarized in Table II. A maximum of 9 partial-waves 
(^Si, ei, ^Po,i.2, ^D2,3, and ^F3.4) were searched in the single-energy analyses. The remaining 
partial-waves were fixed at the energy- dependent values. 

Several Coulomb correction schemes were tested. Our results were found to be relatively 
independent of the chosen form. Differences are most apparent in charge-asymmetry observ- 
ables, as shown in Ref. [^. In the present analysis, we have adopted Hiroshige's formulation 



V. RESULTS AND COMPARISONS 

Results for the partial-wave amplitudes, defined in Eq. (1), are shown in Fig. 3. Over 
our energy range, the dominant amplitudes are the '^Si, '^Pi, '^P2, '^D2, '^Ds, and '^F4. The 
compatibility of our energy- dependent and single-energy solutions is also evident in Fig. 3. 
The lack of single-energy solutions beyond 300 MeV is due to the data distribution displayed 
in Fig. 2. 

As mentioned in Section IV, the energy-dependent solution gives a reasonable overall fit 
to the data. While the x^/datum was about 2 for the selected data base, a much higher 
value would have resulted from the total set of measurements. Some of the data confiicts 
are apparent in Fig. 4, where we have given predictions for observables at Tt^+ = 256 MeV. 
The data are generally well described at this energy, except possibly the most forward iTu 
measurements. We should mention that n'd measurements account for only about a quarter 
of the total dataset. In addition, there is a noticable difference in the x^/datum (1.86 for 
TT^d versus 2.55 for 7c~d). Various n^d total cross sections are plotted in Fig. 5. 

We have also compared our results to those from several other groups. The last single- 
energy partial-wave analyses for vrd elastic scattering were published by the Grenoble- 
Rehovot M, Osaka p[. Saskatoon M, and Mexico-Karlsruhe M groups. These results 
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have generally covered a more narrow energy interval (7 points from 82 to 292 MeV in [Q, 
8 points from 114 to 325 MeV in |[, 9 points from 117 to 324 MeV in 0, and 2 points 
at 256 and 294 MeV in [0). The energy-dependent analysis of Ref. ^ covered the region 
between 65 MeV and 275 MeV. We have plotted our energy-dependent results along with 
several analyses in Fig. 6. The figure displays the experimental version of model (iv) from 
Ref. 0], and the version of Ref. [§ using Blankleider-Afnan |^0| model amplitudes as input. 
The agreement is qualitative at best. 

VI. SUMMARY AND CONCLUSIONS 

In this work we have analyzed a vrd elastic scattering data base which is significantly 
larger than those used in previously published analyses. A reasonable fit to this database 
was found. While there are considerable differences between the results of this and pre- 
vious analyses, one common qualitative feature is present. Several dominant partial- wave 
amplitudes display a "resonancelike" behavior. The correlation between these amplitudes 
is particularly evident in the Argand plot of Fig. 7. This behavior is very similar to that 
found in our recent analyses of it d -^ pp ^ and pp ^ elastic scattering data. 

The present work is completely free of model-based constraints. Previous analyses have 
generally used theory as a guide where insufficient experimental information was available. 
This factor is a likely source for some of the discrepancies evident in Fig. 6. The evolving 
data base provides another. 

While we have analyzed data to 500 MeV, the results above 300 MeV should not be taken 
too seriously. Most of the data base is concentrated below 300 MeV, with only about 200 
measurement covering the region between 300 MeV and 500 MeV. Clearly, we require much 
additional (and more consistent) data to define a unique solution (both above and below 
300 MeV). We should mention that the present solution predicts a rich structure for many 
of the spin-transfer observables defined in Eqs. (3). While only a few of these quantities 
have been measured, some new PSI measurements ||B^ of tjj^, itjo, it2i, and tH at 134, 180, 



and 219 MeV will soon be available. 

This reaction is now incorporated into the SAID program ||6^, which is maintained at 
Virginia Tech. Detailed information regarding the database, partial-wave amplitudes and 
observables may be obtained either interactively, through the SAID system (for those who 
have access to TELNET), or directly from the authors. 
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Figure captions 

Figure 1. Data accumulation from 1952 to the present. Arrows indicate the year when 
measurements of particular observables were first published. 

Figure 2. Energy-angle distribution of total dataset. (a) differential cross section, 
da/dQ for 7r+d, (b) differential cross section, da/dQ for 7r~d, (c) deuteron vector 
analyzing-power, iTu for vr+d, (c) deuteron tensor analyzing-power, T20 for vr+d. 

Figure 3. Partial-wave amplitudes from to 500 MeV. Solid curves are the real parts 
of amplitudes; dashed curves are the imaginary parts. Single-energy solutions 
are plotted as black circles (real part) and open circles (imaginary part). All 
amplitudes have been multiplied by a factor of 10^, and are dimensionless. 

Figure 4. Predictions for vr+d observables at T^^ = 256 MeV. Data have been normal- 



ized, (a) da/dn, (b) zTn, (c) T20, (d) T21, (e) T22, (f) rsi, (g) r22, (h) t 



lab 
20 • 



Figure 5. Total cross sections, a, for n^d scattering, (a) unpolarized total cross section 
ar (solid curve), (b) contribution ^ from n^d — > pp a^ (dashed curve), (c) the 
total elastic cross section erf? (dotted curve), and the remainder Aa given by 

Figure 6. Plots of the amplitudes (a) ^Si, (b) ^Pi, (c) ^¥2, (d) ^T)2, (e) ^Dg, and (d) ^Y^. 
The solid (dashed) curves give the real (imaginary) part from the present energy- 
dependent solution. The circles, squares, triangles, and diamonds are results 
from Refs. 0, 0, [0], and 0, respectively. Filled symbols give real parts; open 
symbols give imaginary parts. All amplitudes have been multiplied by a factor 
of 10'^, and are dimensionless. 

Figure 7. Argand plot of the dominant vrd partial- wave amplitudes '^P2, ^Da, and ^F4 
which correspond to the ^D2, ^Fs, and ^64 pp states, respectively. (Compare 
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Fig. 7 of reference 0). The "X" points denote 50 MeV steps. All amplitudes 
have been multiplied by a factor of 10^. 
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TABLE I. Number and type of observables used in the present analysis to 500 MeV in 
the pion laboratory kinetic energy. The number of excluded data is also given. "^See Refs. 
and |2|. *See Refs. |5[, |g, Q, and |39|. ^See Refs. 0, [|T9[, [gg, and gl[. '^See 
Refs. ig], 13], and [||. ^See Ref. |51|]. ^See Ref. [gg. ^See Refs. |2| and g. 



TT + d TT d 

No. of Data Deleted data No. of Data Deleted data 



Observable 



ax 


57 





69 


2" 


da/dVt 


572 


56^ 


268 


32^ 


a^ 


3 
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3^ 


iTii 


285 


5^ 


6 


1/ 


T20 


42 











T-21 


47 











T22 


76 











Aab 
(-on 


65 


355 









Total 
No. of Energies 



1147 
239 



96 
16 



349 

121 



38 
11 
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Table II. Single-energy (binned) fits and x^ values. Nprm is the number of amplitudes 
(real + imaginary) varied in the fit. Xc ^^ ^^^ ^'^ ^^^ amplitude constraints, Xz) is the 
contribution from data, and x|; is given by the energy-dependent fit, SM94. 



Ttt 

(MeV) 


Range 

(MeV) 


Ndata 
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xl 


xl 


xl 


65 


58.0 - 72.0 


54 


4 


2.4 


86.6 


129.9 


87 


72.0 - 85.5 


24 
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0.4 


20.4 


30.3 


111 


107.5-125.2 


82 
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0.7 


68.9 


68.9 


125 


115.0-134.0 


170 


12 


5.1 


154.1 


205.0 


134 


124.0 - 142.8 


258 


12 


4.2 


293.6 


362.1 


142 


133.0-152.0 


284 


14 


8.1 


345.4 


442.6 


151 


141.0-160.6 
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16 


11.0 


186.9 


256.3 


182 


174.0 - 189.5 
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16 


19.6 


300.3 


445.8 


216 


206.0 - 220.0 
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18 


1.6 


121.3 


148.1 


230 


220.0 - 238.0 


53 


18 


7.8 


51.3 


96.5 


256 


254.0 - 260.0 
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18 


2.9 


132.5 


200.6 


275 


270.5 - 284.4 


40 


18 


2.9 


17.4 


42.6 


294 


284.4 - 300.0 


132 


16 


3.3 


228.5 


263.6 
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